ARBRL-TR-2031 


Ab-AK  I 6S| 


AD 


TECHNICAL  REPORT  ARBRL-TR-2031 

AIR  DEFENSE  GUN  FIRE  CONTROL  SYSTEM 
PERFORMANCE  WITH  A BIAS  PREDICTOR 

PakT.  Yip 

December  1977 


Approved  for  public  release;  distribution  unlimited. 


ffUS  QtrAlIW  niSPEWEK  B 


USA  ARMAMENT  RESEARCH  AND  DEVELOPMENT  COMMAND 
USA  BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GROUND,  MARYLAND 


Destroy  this  report  when  it  is.no  longer  needed. 
Do  not  return  it  to  the  originator. 


Secondary  distribution  of  this  report  by  originating 
or  sponsoring  activity  is  prohibited. 

Additional  copies  of  this  report  may  be  obtained 
from  the  National  Technical  Information  Service, 

U.S.  Department  of  Commerce,  Springfield,  Virginia 

22161. 


The  findings  in  this  report  are  not  to  be  construed  as 
an  official  Department  of  the  Army  position,  unless 
so  designated  by  other  authorized  documents. 


The  use  jf  tmde  ruvrtes  ov  manufacturevs ' nameo  in  thin  r*eport 
does  not  oonatitute  indorsement  of  any  aormeraial  -product. 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entarad) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER 

TECHNICAL  REPORT  ARBRL-TR-2031 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  {'and  Subtif/e) 

AIR  DEFENSE  GUN  FIRE  CONTROL  SYSTEM 
PERFORMANCE  WITH  A BIAS  PREDICTOR 

5.  type  of  REPORT  & PERIOD  COVERED 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7,  AUTHORr*; 

PAK  T.  YIP 

8.  CONTRACT  OR  GRANT  NUMBERfs; 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

US  Army  Ballistic  Research  Laboratory 

Aberdeen  Proving  Ground,  MD  21005 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  4 WORK  UNIT  NUMBERS 

1L662618AH80 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  Materiel  Development  ^ Readiness  Command 
5001  Eisenhower  Avenue 

Alexandria,  VA  22333 

12,  REPORT  DATE 

DECEMBER  1977 

13.  NUMBER  OF  PAGES 

34 

14.  MONITORING  AGENCY  NAME  & ADD RESS(if  different  from  Controlling  Office) 

15.  SECURITY  CLASS,  (of  thia  report) 

Unclassified 

15a,  DECLASSIFIC  ATI  on/ DOWNGRADING 

SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  thia  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abatract  entered  in  Block  20,  if  different  from  Report) 


\ 


18.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  reverse  aide  if  necessary  and  identify  by  btock  number) 


Gun  Air  Defense  Fire  Control  System  FACTS 

Kalman  filter  Bias  Predictor 

Probability  of  hit 


20.  ABSTRACT  f Continue  on  reverse  aide  if  neceaaaay  and  identify  by  block  number) 

A bias  predictor  was  introduced  to  couple  with  a simple  nine  state  discveti^ 
Kalman  filter  for  solving  the  gun  fire  control  problem.  The  performance  of 
this  system  showed  that  an  earlier  and  higher  probability  of  hit  can  be 
achieved  if  the  defended  target  position  is  known. 


DD  1 1473  EOrTlOH  OF  T NOV  6S  IS  OBSOLETE 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


TABLE  OF  CONTENTS 


Page 

LIST  OF  ILLUSTRATIONS 5 

I.  INTRODUCTION 7 

II.  THE  FILTER  MODEL 9 

III.  THE  BIAS  PREDICTOR H 

IV.  THE  SYSTEM  PERFORMANCE  MODEL 14 

V.-  RESULTS 19 

VI.  DISCUSSION 30 

REFERENCES 31 

DISTRIBUTION  LIST 33 


3 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

1 Typical  Estimation  Error  Profile  for  Tracking  an 

Attacking  Aircraft  8 

* 2 Point  Mass  Force  System  for  an  Aircraft 13 

3 An  Aircraft  Attacking  Flight  Profile  18 

4 Projectile-Aircraft  Miss  Distance  Against  Aircraft 

Range  at  Intercept  for  FACTS  Profile  as  Shown  in 

Figure  3 19 

5 Probability  of  Hit  per  Burst  (33  Rounds]  With  Gun  at 

Origin  in  FACTS  Profile  as  Shown  in  Figure  3 21 

6 Probability  of  Hit  per  Burst  (33  Rounds)  With  Gun 

Displaced  by  1 km  Cross  Range  in  FACTS  Profile  as 

Shown  in  Figure  3 22 

7 Probability  of  Hit  per  Burst  (33  Rounds)  With  Gun 

Displaced  by  2 km  Up  Range  in  FACTS  Profile  as 

Shown  in  Figure  3 23 

8 Aircraft  Attacking  Profile  24 

9 Probability  of  Hit  per  Burst  (33  Rounds)  With  Gun  at 

Origin  in  FACTS  Profile  as  Shown  in  Figure  8 25 

10  Probability  of  Hit  per  Burst  (33  Rounds)  With  Gun 

Displaced  by  2 km  Up  Range  in  FACTS  Profile  as 

Shown  in  Figure  8 26 

11  Aircraft  Attacking  Profile  27 

12  Probability  of  Hit  per  Burst  (33  Rounds)  With  Gun  at 

Origin  in  FACTS  Profile  as  Shown  in  Figure  11  ....  28 

13  Probability  of  Hit  per  Burst  (33  Rounds)  With  Gun  at 

“■  Displaced  by  2 km  Up  Range  in  FACTS  Profile  as 

Shown  in  Figure  11 29 


5 


I . INTRODUCTION 
' ■ ■ ■ ■' 

In  general,  the  strategy  of  an  attacking  aircraft  may  consist  of 
two  major  typical  phases,  the  pre-delivery-maneuver  phase  followed  by 
the  weapon-delivery  phase.  In  the  pre-delivery-tnaneuver  phase,  the  air- 
craft usually  maneuvers  evasively  creating  a highly  jinking  noise  signal 
with  which  very  poor  estimation  of  the  aircraft  states  can  be  expected. 

A typical  estimation  error  profile  is  demonstrated  in  Figure  1.  It  is 
obvious  that  the  poor  state  estimates,  resulting  from  this  maneuver  phase, 
will  present  poor  intercept  prediction  of  either  second  order  or  first 
order  type.  However,  a predictor  making  use  of  information  on  aircraft 
position,  velocity  and  ground  target  location,  and  eliminating  the  esti- 
mates of  aircraft  acceleration  and  jerk  might,  in  concept,  provide  more 
accurate  prediction  capability.  That  is,  a bias  predictor  making  use 
mainly  of  the  position  estimates  and  estimated  speed  of  the  aircraft,  and 
the  assumption  of  a known  target  position,  may  provide  a better  intercept 
prediction.  In  the  weapon-delivery  phase,  the  aircraft  uses  only  those 
minimal  maneuvers  required  to  maintain  target  alignment.  An  earlier  hit 
in  this  phase  has  normally  been  limited,  due  to  the  slow  settling  of  the  ve- 
locity and  acceleration,  state  estimates  of  the  aitcraft,  changing  from  a 
highly  jinking  to  a slight  maneuver  phase.  We  will  attempt  to  use  only 
the  position  and  speed  estimates  of  the  aircraft  and  knowledge  of  target 
position  to  generate  anl  intercept  predictidn,  hoping  that  will  help  achieve 
earlier  hits.  Also,  it  is  in  this  phase  that  we  may  use  an  appropriate 
model  with  good  aerodynamic  characteristics  to  simulate  the  attacking 
aircraft. 

Therefore,  this  report  is  devoted  to  the  study  of  the  possibility 
of  improving  system  performance  by  using  the  bias  predictor  for  intercept 
prediction  in  the  weapon-delivery  phase.  The  extension  of  the  idea  to 
the  aircraft  prediction  in  the  pre-delivery  phase  will  be  examined  in  a 
follow-up  study.  ■; 

In  Section  2,  a Kalman  filter  with  an  adaptive  scheme  for  the  plant 
noise  covariance  matrix  is  discussed.  A bias  predictor  using  aircraft 
aerodynamics  is  presented  in  Section  3.  In  Section  4,  a system  perform- 
ance model  is  discussed.  Finally,  some  results  from  applying  the  Kalman 
filter  and  the  bias  predictor  on  several  passes  of  the  FACTS*  profile  are 
summarized  in  Section  5. 


*”Frankford  Aircraft  Capability  Test  Study, actual  flight  data  prepared 
by  Frankford  Arsenal,  Department  of  the  Army,  Philadelphia,  PA,  March, 
1974.  . , 
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Figure  1.  Typical  Estimation  Error  Profile  for  Tracking  an  Attacking  Aircraft 


II.  THE  FILTER  MODEL 


A simple  9 state  discrete  Kalman  filter  was  used  to  provide  estimated 
aircraft  states.  The  plant  equation  describing  aircraft  dynamics  is: 

t 

Y = T Y + T1 

K+1  'k+1,K  K K 


where  X is  the  state  vector  at  time  KAt;  X, 


the  state  vector  at  time 


(K+1) the  state  transition  matrixf^i^^,  the  driving  input;  and 
At,  the  sampling  time  step.  The  estimated  state  equation  is 


X- 


K+1 


/\ 

T Y + W 

K+1,K  K,  K+1 


K+1 


"k.i  ’’k.i 


,K 


where  X^^  is  the  estimated  state,  vector  at  time  KAt;  ^ the  estimated 

state  vector  at  time, (K+1) At ; the  measurements  ar  time  (K+1) At; 

Hk+i^  the  observation  matrix;  and  Kalman  gain  matrix. 

The  gain  matrix,  W,  is  a result  of  state  error  propagation  which  is 
a function  of  the  plant  noise  covariance  matrix  defined  as 


Q = Expected  Value  of 

K > ■ 


Following  Singer's  derivation,  this  plant  noise  covariance  has  a 
form  when  At  is  sufficiently  small  so  that  aAt<<%, 


Q = 2aA  t (7^ 
^ m 


2.  9-  — 

£ 0 o_ 

o 0 I 


where  « is  related  to  the  maneuver  time  constant;  <r^  , the  variance  of 
aircraft  acceleration;  o_,  a 3 by  3 zero  matrix;  and  a 3 by  3 identity 
matrix. 

2 


And, 


where  maximum  acceleration  the  aircraft  may  use;  Fp  ^ 

a function  of  and  P^;  Pj^^  the  probability  that  the  aircraft  is  flying 


^Robert  A.  Singer,  ^’Evaluating  Optimal  Tracking  Filter  Performance  for 
Manned  Maneuvering  Targets,”  IEEE  Trans . Aerospace  and  Electronic  Systems, 
Vol.  AES-6,  No.  4,  July  1970. 
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with  an  acceleration  probability  that  the  aircraft  is 

flying  with  constant  velocity. 

It  is  well  known  that  this  plant  noise  covariance  is  strongly  related 
to  estimation  accuracy.  For  a given  maneuvering  aircraft  and  a fixed  pro- 
cessing data  rate,  values  of  Qj(  that  are  too  large  will  pass  excessive 
measurement  noise,  and  values  of  Q](  too  small  will  block  off  the  critical 
information  of  aircraft  maneuvers.  Some  designers  advocate  bandwidth 
adaptation.^  For  this  study,  the  guesses  of  acceleration  probability  are 
adaptive;  hence,  the  variance  is  adaptive. 

To  accomplish  this  scheme,  five  intervals  of  radial  acceleration  rate 
are  predefined.  Five  sets  of  acceleration  probabilities  corresponding  to 
these  intervals  are  assigned.  The  average  of  the  previous  ten  estimates 
of  aircraft  acceleration  is  stored.  The  rate  of  change  of  this  average 
acceleration  is  tested  against  the  intervals  and  the  acceleration  proba- 
bilities are  updated  accordingly.  This  filter,  similar  to  a hypothesis 
testing  or  an  adaptive  bandwidth  filter,  has  been  studied  and  a report  is 
being  prepared  to  be  published  soon.  The  study  of  this  filter  has  indi- 
cated an  improvement  of  about  13  percent  in  estimation  accuracy  over  that 
with  a fixed  value  plant  noise  covariance. 

As  aircraft  state  estimates  are  obtained  from  this  filter,  the  first 
and  second  order  predictions  are  computed  by  the  following  equations  re- 
spectively: 


A A 
^IST  ^ ^ 


_ A A , A 2 

A A , . , , . 

where  X is  the  estimated  position  vector;  X^,  the  estimated  velocity 

vector;  the  estimated  acceleration  vector;  and  tp,  the  projectile 

time  of  flight  to  intercept  which  is  computed  using  the  estimated  states 

and  a projectile  flight  equation  in  a performance  model. 


2 

C.  M.  Brown,  Jr.,  and  C.  F.  Price,  "Adaptive  Tracking  Filter  Design  and 
Evaluation  for  Gun  Fire  Control  Systems,"  Reading,  MA,  The  Analytical 
Science  Corporation,  January  1974. 
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III.  THE  BIAS  PREDICTOR 


In  building  a bias  predictor,  we  limited  our  interest  to  the  fol- 
lowing conditions.  The  target  position  is  known  to  the  fire  control. 

The  aircraft  will  attack  the  target  with  a weapon-delivery  phase.  A 
nine  state  Kalman  filter  will  provide  the  aircraft  position,  velocity, 
and  acceleration  estimates  as  inputs  to  this  predictor. 

The  bias  predictor  performs  the  following:  (1)  process  the  inputs 
to  identify  aircraft  strategy,  (2)  set  initial  and  terminal  conditions 
for  flying  a simulated  aircraft,  and  (3)  integrate  the  simulated  aircraft 
dynamic  equations  through  projectile  time  of  flight. 

- Processing  the  inputs  is  no  more  than  computing  the  bearing  of  the 
aircraft  relative  to  its  target,  and  the  initial  pitch  angle  of  the  air- 
craft. Assuming  the  aerodynamic  angle  of  attack  to  be  small,  with  the 
information  of  target  position,  aircraft  position  and  velocity,  the  above 
computation  is  just  simple  vector  geometry. 


To  form  the  rules  of  the  game,  an  aircraft  roll  is  called  for  when- 
ever the  magnitude  of  the  relative  bearing  is  larger  than  2 , and  a climb 
or  dive  is  demanded  whenever  the  pitch  angle  is  greater  than  3 or  less 
than  -3  respectively.  According  to  these  rules,  types  of  aircraft  strat- 
egy are  identified  as  positive  pitch,  positive  pitch  and  roll,  horizontal 
flight,  horizontal  flight  and  roll,  dive,  and  dive  and  roll.  There  is  an 
additional  type  which  is  an  approximate  180  roll  and  dive  switched  from 
the  positive  pitch  type  when  the  pitch  angle  is  greater  than  45  . 


The  initial  conditions  are  set  with  the  inputs  and  as  functions  of 
type  of  aircraft  strategy.  The  initial  estimated  position  and  velocity, 
the  pitch  angle,  and  the  heading  angle  of  the  aircraft  are  given  as  the 
inputs  or  computed  directly  from  the  inputs.  The  acceleration,  roll 
angle,  and  percentage  of  power  in  use  are  initialized  as  functions  of 
type  of  aircraft  strategy.  Maximum  values  for  acceleration,  roll  angle, 
roll  rate,  rate  of  acceleration,  and  rate  of  power  are  fixed  as  functions 
of  type  of  aircraft  strategy  also.  The  initial  conditions  coincide  with 
the  initial  estimation  of  the  aircraft  state.  However,  the  terminal 
conditions  demand  the  simulated  aircraft  to  fly  a certain  bias  profile. 


The  terminal  conditions  are  the  following:  The  commanded  pitch 
angle  is  a function  of  aircraft  strategy  and  the  initial  pitch  angle. 
The  attack  altitude  is  a function  of  aircraft  strategy.  The  commanded 
velocity  is  a function  of  the  commanded  pitch  angle.  In  the  case  of 
the  attack  heading  for  the  simulated  aircraft,  a parametric  formula  was 
used. 
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Attack  heading  = Initial  aircraft  heading  + (1  + ) * Relative  bear- 

ing of  aircraft  and  its  targets,  o 

where  K is  an  arbitrary  constant,  V is  the  estimated  speed  of  the  air- 
craft, and  V is  set  to  Mach  one.  The  constant  K is  generally  set  to 
0.3.  ° 

3 

An  Air  Force  Model,  the  Blue  Max  , was  used  to  simulate  aircraft 
aerodynamics.  Basically,  a set  of  general  point  mass  force  equations 
is  used  in  a body  fixed  right-handed  coordinate  system.  The  dynamic 
equations  are: 


X = (T-D-Mg  sin  0)  /M 

y = N sin  '^/M 

*z  = (N  cos  0 - Mg  cos  0)/M 

where  x is  the  aircraft  acceleration  along  the  body  axial  axis;  y,  that 
along  the  horizontal  axis;  z,  that  along  the  positive  upward  axis;  T,  the 
thrust  exerted  on  the  aircraft;  D,  the  axial  drag-,  N,  the  normal  force; 

M,  the  mass  of  the  aircraft;  g,  the  gravitation  force  constant;  0,  the 
pitch  angle;  and  0,  the  roll  angle.  The  system  is  shown  in  Figure  2. 

The  driving  rate  controls  of  the  aircraft  have  the  following  forms: 

p = (h^  - h)  - h 

V - V)  - V 
G = (0^  - 0)  - q 

where  p is  the  angular  roll  rate;  h,^the  present  heading;  h^,‘the  com- 
manded heading;  h,  the  heading  rate;  T^,  the  throttle  rate;  V,  the^pre- 
sent  speed;  V , the  commanded  speed;  the  rate  change  of  speed;  G,  the 
G-load  rate;  6,  the  pitch  angle;  0^,  the  commanded  pitch  angle;  and  q, 
the  pitch  rate. 

A power  table  is  provided  as  a function  of  altitude  and  speed  of  the 
aircraft  for  the  computation  of  thrust.  A table-look-up  procedure  is 
used  for  computing  the  normal  force  as  a function  of  aircraft  speed.  A 
drag  table  is  provided  as  a function  of  altitude  and  speed  of  the  air- 
craft for  obtaining  the  drag  force.  With  these  forces  and  computed  pitch 


3 

Oleg  Komarnitsky  and  AF/SAFG  personnel,.  ”Blue  Max,’*  program  user’s  manual, 
ACS/Studies  and  Analysis  Office,  Headquarters  USAF,  1974. 
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Figure  2.  Point  Mass  Force  System  for  an  Aircraft 
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angle,  roll  angle,  and  the  aircraft  heading,  the  dynamic  equations  are 
integrated  through  one  time  step.  The  process  is  repeated  through  the 
given  projectile  time  of  flight.  Then,  the  predicted  position  of  the 
aircraft  from  this  bias  predictor  is  the  actual  future  position  of  a sim- 
ulated aircraft. 


IV.  THE  PERFORMANCE  MODEL 


To  evaluate  the  performance  of  this  gun  fire  control  system,  a 
modified  STAG"^  model  was  used.  Basically,  the  STAG  computes  the  time  of 
flight  of  the  projectile  to  intercept,  the  mean  distance  at  closest  ap- 
proach, and  the  probability  of  hit  against  an  aircraft  with  a burst  of 
n projectiles.  The  STAG  model  provides  an  analytic  solution  with  a sin- 
gle run,  but  it  is  not  designed  to  include  second  order  prediction  or 
system  lag.  Bias  and  dispersion  are  taken  only  root  sum  square  terms 
from  input.  To  modify  the  STAG,  the  probabilities  of  events  were  rede- 
fined. A new  error  computation  routine  was  added. 

The  probability  of  hit  against  a diffused  target  with  a burst  of 
n rounds  is  given  as  in  the  STAG  model,  by 


j-1 


V 


where  Ay  is  the  projected  vulnerable  area  of  the  target;  h,  the  mean 
distance  at  closest  approach;  <r  , the  standard  deviation  of  "random" 
bias;  and  <r^,  the  standard  deviation  of  round  to  round  dispersion. 

In  general,  the  miss  distance  is  a function  of  the  estimated  states 
and  the  prediction  algorithm.  The  errors  embedded  in  the  estimated  states 
are  system  lag,  sensor  random  noise,  aircraft  glint,  and  modeling  noise 
or  maneuver  noise. 

In  consideration  of  the  randomness  of  the  measurements  for  a partic- 
ular flight  profile,  the  average  P for  20  runs  was  used  as  a system  per- 
formance index.  The  starting  values  of  the  random  number  generator  for 


"^Fred  Bunn,  "Surface  to  Air  Gun  Model,"  unpublished  study  report.  Air 
Warfare  Division,  AMSAA,  APG,  MD,  1971. 
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sensor  noise  and  aircraft  glint  were  different  run  by  run.  The  error  bud- 
get used  for  this  performance  analysis  was  based  on  information  from  con- 


tractors as  follows: 

Sensor  noise,  azimuth  (la) 00075  radians 

elevation  00075  radians 

I'ange 3 meters 

Aircraft  glint,  azimuth 1 meter 

elevation 1 meter 

Wind,  northern 5 m/s 

eastern 5 M/s 

Gun  recoil,  azimuth 0005  radians 

elevation 0005  radians 

Gun  boresight,  azimuth 001  radians 

elevation 001  radians 

Sensor  misalignment,  azimuth 00075  radians 

elevation 00075  radians 

Sensor  range  bias ' . . . . 2 meters 

Muzzle  velocity  variation 6 M/s 

Projectile  dispersion,  transverse.  . .001  radians 

elevation  . . .001  radians 


An  additional  factor  in  the  error  budget  was  the  artificial  disper- 
sion which  is  mostly  a function  of  range  and  aircraft  acceleration.  It 
may  be  in  the  form 


where  cr  is  the  standard  deviation  of  artificial  dispersion;  cr  the 

nonlinear  range  dependent  portion;  C,  a constant;  and  X , the  estimated 
acceleration  vector.  ^ 
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If  the  statistics  of  aircraft  maneuvers  are  completely  known  to  the 
gun  fire  control,  artificial  dispersion  is  not  desirable  for  an  optimal 
predictor.  In  the  case  that  neither  the  statistics  of  the  attacking 
aircraft  maneuvers  are  completely  known,  nor  the  distribution  of  the  air- 
craft acceleration  is  exactly  Gaussian  as  modeled,  appropriate  artificial 
dispersion  often  helps  improve  the  system  performance.  Since  the  purpose 
of  this  study  was  to  examine  the  possibility  of  improving  system  perform- 
ance by  using  a bias  predictor,  optimized  artificial  dispersion  was  used. 
In  this  context,  all  predictors  in  this  paper  were  compared  in  the  upper 
bound  of  system  performance.  It  has  been  found  that  the  optimized  dis- 
persion is  approximately  seven-tenths  of  the  miss  distance  at  closest 
approach. 


With  regard  to  contribution  to  miss  from  various  error  sources, 
readers  are  referred  to  the  report  by  H.  K.  Weiss. ^ In  his  report, 
angular  errors  due  to  wind,  time  of  flight  variation,  and  slant  range 
bias  are  derived.  Muzzle  velocity  variation  and  head  wind  give  rise  to 
time  of  flight  variation  which  in  turn  contributes  to  angular  miss. 


With  the  filter,  the  predictors,  and  the  performance  model,  we  need 
a projectile  and  some  flight  profiles  to  play  the  gun  air  defense  game. 

The  following  values  from  a typical  modern  gun  system  are  assumed. 
Muzzle  velocity  was  set  to  1176  meters  per  second.  The  drag  factor 
^ was  set  to  0.129.  Then,  the  projectile  time  of  flight 

t = R 

F V - /3R 

m ^ 

where  R is  the  intercept  range.  Thirty- three  rounds  per  burst  was 
chosen  and  the  time  between  burst  was  2 seconds.  Some  flight  profiles 
were  selected  from  the  FACTS  results.  The  FACTS  data  were  then  processed 
by  least  square  polynomial  curve  fitting.  Every  eleven  points  were  pro- 
cessed to  obtain  five  with  a view  to  maintaining  the  moderately  high  fre- 
quency component  of  aircraft  motion.  The  smoothed  data  were  then  cor- 
rupted by  Gaussian  noise  and  fed  into  the  filter. 


^Harry  L.  Reed,  Jr.,  "Some  Bounds  on  the  Generalized  Fire  Control  Problem," 
Report  No.  1946,  Ballistic  Research  Laboratory,  APG,  MD,  November,  1976. 
f:AD  #A033043 

H.  K.  Weiss,  "Closed  Loop  Systems  Concept  Study,"  final  report.  Data 
Systems  Division,  Litton  Systems,  Inc.,  Van  Nuys,  CA,  September,  1974. 
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V.  RESULTS 


For  the  FACTS  flight  profile,  as  shown  in  Figure  3,  miss  distance 
versus  target  range  at  intercept  for  the  first  order,  second  order,  and 
bias  predictors  is  shown  in  Figure  4.  The  results  are  summarized  in 
Table  1 and  shown  in  Figure  5. 

With  the  same  flight  path,  except  that  the  fire  control  was  placed 
1000  meters  cross  range  in  one  case,  and  2000  meters  up  range  in  another 
case,  the  results  are  shown  in  Figures  6 and  7 respectively. 

Results  for  the  profile  in  Figure  8 are  shown  in  Figures  9 and  10. 
For  the  profile  in  Figure  11,  results  are  shown  in  Figures  12  and  13. 


Figure  3.  An  Aircraft  Attacking  Flight  Profile 
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Figure  4.  Projectile-Aircraft  Miss  Distance  Against  Aircraft  Range 
at  Intercept  for  FACTS  Profile  as  Shown  in  Figure  3, 


TABLE  1.  SUMMARIZED  PERFORMANCE  OF  SYSTEMS  WITH  DIFFERENT 

PREDICTORS  FOR  FLIGHT  PROFILE  AS  SHOWN  IN  FIGURE  3 


POINT  OF 
INTERCEPT 

1 

2 

3 

4 

5 

r 

6 

7 

AVERAGE 

APPROXIMATE 

I NT.  RANGE (m) 

3535 

3257 

2962 

2646 

2307 

1949 

1571 

APPROXIMATE 

TIME  OF  FLIGHT 
(s) 

4.909 

4.310 

3.732 

3.170 

2.626 

2.108 

1.614 

MISS  DISTANCE 
(m) 

1st  Order 

2nd  Order 

Bias  Pred. 

321.4 

142.3 

309.8 

82.9 

178.5 

40.7 

68.0 

170.9 

52.8 

11.7 

30.5 

28.6 

9.8 

14.4 

7.3 

5.5 

8.0 

7.8 

19.2 

16.5 

15.7 

PROBABILITY 

OF  SINGLE 

SHOT  HIT 

1st  Order 

2nd  Order 

Bias  Pred. 

.0000 

.0001 

.0000 

.0001 

.0000 

.0005 

.0001 

.0000 

.0002 

.0050 

.0012 

.0012 

.0094 

.0067 

.0159 

.0208 

.0139 

.0187 

.0030 

.0043 

.0050 

PROBABILITY 

OF  BURST  HIT 
(33  Rounds) 

1st  Order 

2nd  Order 
Bias  Pred. 

.0003 

.0033 

.0002 

,0026 

.0005 

.0155 

.0031 

.0005 

.0056 

.1512 

.0399 

.0393 

.2677 

.1992 

.4102 

.5005 

.3709 

.4628 

.0933 

.1340 

.1526 

i 

.146 

.107 

.155 
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Figure  5.  Probability  of  Hit  per  Burst  (p5  Rounds)  With,  Gun  at  Origin 
in  FACTS  Profile  as  Shown  in  Figure  3. 
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Figure  6.  Probability  of  Hit  per  Burst  (33  Rounds)  With  Gun  Displaced 
by  1 km  Cross  Range  in  FACTS  Profile  as  Shown  in  Figure  3. 
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Figure  7.  Probability  of  Hit  per  Burst  (33  Rounds)  with  Gun  Displaced 
by  2 km  Up  Range  in  FACTS  Profile  as  Shown  in  Figure  3. 
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Pigui'e  9.  Probability  of  Hit  per  Burst  (33  Rounds)  With  Gun  at  Origin 
in  FACTS  Profile  as  Shown  in  Figure  8, 
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Figure  10.  Probability  of  Hit  per  Burst  (33  Rounds]  With.  Gun  Displaced 
by  2 km  Up  Range  in  FACTS  Profile  as  Shown  in  Figure  8, 
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Figure  11.  Aircraft  Attacking  Profile. 
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Pigu;]70  12.  Pirobs-bility  of  Hit  pox  Biixst  (33  Rounds}  Witb  Gun  nt  Oxigin 
in  FACTS  Pxofile  as  Shown  in  Figuxe  11. 
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Figure  13.  Probability  of  Hit  per  Burst  (33  Rounds)  With  Gun  Displaced 
by  2 km  Up  Range  in  FACTS  Profile  as  Shown  in  Figure  11. 


VI.  DISCUSSION 


In  general,  a Kalman  filter  is  always  included  in  the  design  of 
modern  gun  fire  control  systems.  Since  an  attacking  aircraft  has  a large 
spectrum  of  possible  acceleration  rates,  the  Kalman  prediction  has  cer- 
tain disadvantages.  In  the  case  of  the  second  order  prediction,  very 
little  accuracy  can  be  obtained  in  the  beginning  stage  of  the  weapon- 
delivery  phase  because  the  filter  requires  approximately  3 seconds  to 
settle  its  acceleration  states  from  high  maneuver  to  very  low  maneuver 
at  ranges  beyond  3 kilometers.  Also,  for  the  relative  slight  maneuvers 
in  this  phase,  the  processing  of  acceleration  noise  by  the  filter  will  re- 
sult in  an  increase  in  the  prediction  error.  Therefore,  in  this  case,  the 
second  order  prediction  will  be  inferior, to  the  first  order  prediction. 

In  the  case  of  the  first  order  prediction,  about  one  second  is  required 
for  the  velocity  term  to  settle  for  the  same  transition  of  maneuver. 
Moreover,  if  the  flight  profile  is  a dive  with  an  attack  elevation  angle, 
prediction  errors  may  be  increased  due  to  neglecting  the  gravitational 
pull  on  the  aircraft. 

In  this  study,  it  has  been  demonstrated  that,  if  the  aircraft  act- 
ually attacks  a target,  the  position  of  which  is  known,  a bias  prediction 
algorithm  provides  better  state  predictions  in  the  weapon-delivery  phase. 
This  predictor  becomes  effective  as  soon  as  reasonable  position  estimates 
are  obtained.  The  bias  predictor  model  describes  the  aircraft  profile 
more  accurately  in  lay  down  or  delivery  dive  because  the  flight  is  aero- 
dynamically  simulated.  However,  this  bias  predictor  adds  a considerable 
burden  to  the  fire  control  computer. 

As  now  there  are  still  great  difficulties  in  providing  accurate 
estimation  of  aircraft  states  in  the  pre-delivery-maneuver  phase  in  the 
attacking  aircraft  strategy.  No  known  air  defense  fire  control  algorithm 
is  adequate  for  this  task,  yet,  as  aircraft  fire  control  systems  become 
more  sophisticated,  their  capability  for  delivery  of  ordnance  at  greater 
ranges  and  maneuver  conditions  will  also  increase.  It  is  this  unresolved 
problem  that  creates  the  challenge  for  further  research  effort.  There- 
fore, our  next  study  will  be  addressed  to  using  an  aircraft  model  with 
identified  sets  of  aerodynamic  characteristics  of  corresponding  aircraft 
to  predict  projectile-aircraft  intercept  in  the  pre-delivery -maneuver 
phase. 

To  predict  the  future  position  of  an  aircraft  is  never  an  easy  task, 
since  behind  the  aircraft  a human  decision  is  involved.  Improvement  of 
system  performance  has  come  from  advanced  sensors,  high  velocity  pro- 
jectiles, better  filters  and  better  gun  systems.  Perhaps  switching  of 
predictors  for  the  different  phases  of  an  attacking  flight  profile  may 
join  the  march  on  the  road  with  the  rest  to  a better  solution  for  the 
gun  air  defense  problem. 
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